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ABSTRACT 

Cosmic structures leave an imprint on the microwave background radiation through 
the integrated Sachs-Wolfe effect. We construct a template map of the linear signal 
using the SDSS-III Baryon Acoustic Oscillation Survey at redshift 0.43 < 2 ; < 0.65. 
We verify the imprint of this map on the Planck CMB temperature map at the 97% 
confidence level and show consistency with the density-temperature cross-correlation 
measurement. Using this ISW reconstruction as a template we investigate the presence 
of ISW sources and further examine the properties of the Granett-Neyrinck-Szapudi 
supervoid and supercluster catalogue. We characterise the three-dimensional density 
profiles of these structures for the first time and demonstrate that they are significant 
structures. Model fits demonstrate that the supervoids are elongated along the line-of- 
sight and we suggest that this special orientation may be picked out by the void-finding 
algorithm in photometric redshift space. We measure the mean temperature profiles in 
Planck maps from public void and cluster catalogues. In an attempt to maximise the 
stacked ISW signal we construct a new catalogue of super-structures based upon local 
peaks and troughs of the gravitational potential. However, we do not find a significant 
correlation between these structures and the CMB temperature. 

Key words: cosmological parameters, observations, large-scale structure of Universe, 
methods: statistical 


1 INTRODUCTION 

Striking correlations have been observed between cosmic 
super-structures and the cosmic microwave background 
(CMB). Most prominently Granett, Neyrinck & Szapudi 
(2008) (hereafter GraNeSz) reported a 4cr detection of voids 
and clusters at redshift z = 0.5 identified in the Sloan Dig¬ 
ital Sky Survey (SDSS) imprinted on the CMB as lO/xK 
cold and hot spots on 4° scales. Additionally, there is evi¬ 
dence for a supervoid at redshift z = 0.2 aligned with the 
Cold Spot suggesting a causal relationship (Szapudi et al. 
2015). While these hndings may be statistically significant, 
the physical mechanism producing the correlation has yet to 
be determined. The leading candidate has been the late-time 
integrated Sachs-Wolfe effect (Sachs & Wolfe 1967; Rees & 
Sciama 1968). However, the scale of the temperature corre¬ 
lation is inconsistent with the expectation in the standard 
ACDM cosmological model for typical supervoid and super¬ 
clusters (Cai et al. 2010; Nadathur, Hotchkiss & Sarkar 2012; 
Hernandez-Monteagudo & Smith 2013; Flender, Hotchkiss 
& Nadathur 2013; Nadathur & Hotchkiss 2014; Hotchkiss 
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et al. 2015; Aiola, Kosowsky & Wang 2015) and cannot fully 
explain the Cold Spot (Nadathur et al. 2014; Zibin 2014). 
Moreover, the signal-to-noise ratio for the ISW measurement 
with SDSS galaxies is expected to be ~ 1.7 (Hernandez- 
Monteagudo et al. 2014), while a full-sky survey to redshift 
2 can reach a maximum S/N 7 (Afshordi 2004). Thus the 
anomalous contribution from individual sources may point 
to new physics. 

The ISW signal is challenging to measure because the 
amplitude is an order of magnitude below the primary 
temperature anisotropies. It may be detected by cross- 
correlating galaxy tracers with the cosmic microwave back¬ 
ground over the sky or at the locations of particular super¬ 
structures (Nishizawa (2014) offers a review). In flat cos¬ 
mologies the ISW signal is sensitive to dark energy (Crit¬ 
tenden & Turok 1996) and provides important constraints 
for modified gravity models (Giannantonio et al. 2010) as 
well as for primordial non-Gaussianity (Afshordi & Tolley 
2008). The signal may be enhanced by adding constraints 
from the primordial polarisation of the CMB (Frommert & 
EnBlin 2009). It may also be independently estimated by 
correlating against the infrared background (Ilic et al. 2011) 
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and 21cm brightness maps (Raccanelli et al. 2015) adding 
to the detection significance. 

Since the ISW effect is subtle, analyses can be prone 
to over-fitting and suffer from the ‘look elsewhere’ effect 
(Peiris 2014). It is a particular problem for the study of ISW 
sources since the physical nature of supervoid and superclus¬ 
ter structures identified in galaxy surveys is not well-known 
and there is not a clear expectation for the signal. Numeri¬ 
cal N-body simulations have been used to predict the ISW 
anisotropy arising from voids and clusters and to determine 
the optimal detection filters (Cai et al. 2010, 2014; Hotchkiss 
et al. 2015). Guided by these results, attempts to repro¬ 
duce the GraNeSz measurement with other super-structure 
catalogues have lead to only marginal or null detections 
(Kovacs & Granett 2015; Hotchkiss et al. 2015; Nadathur 
& Hotchkiss 2014; Gai et al. 2014; Ilic, Langer & Douspis 
2013). These studies would suggest that the GraNeSz de¬ 
tection arises from a different physical mechanism or is a 
statistical fluke. 

However, an important piece of the puzzle remains out 
of place. To date, there has been no characterisation of 
the three-dimensional properties of the GraNeSz structures. 
Analyses have used superstructure catalogues derived from 
spectroscopic redshift catalogues employing variations of the 
watershed algorithm based on the ZOBOV code (Neyrinck 
2008; Neyrinck, Gnedin & Hamilton 2005). The GraNeSz 
super-structures were identified with a photometric redshift 
catalogue and the correspondence between the void find¬ 
ing techniques is not yet clear. As pointed out by Kovacs 
V Granett (2015) the GraNeSz supervoids do not correlate 
with the voids identified in three-dimensions with the wa¬ 
tershed algorithm. 

Besides the use of N-body simulations, an alternative 
empirical approach is possible: we may predict the ISW sig¬ 
nal directly from galaxy surveys, without modelling indi¬ 
vidual structures. All-sky reconstructions of the ISW signal 
assuming linear dynamics have already proven to be use¬ 
ful for investigating possible large-scale anomalies on the 
GMB (Maturi et al. 2007; Granett, Neyrinck & Szapudi 
2009; Francis & Peacock 2010; Rassat, Starck & Dupe 2013; 
Manzotti V Dodelson 2014). 

Previously, the GraNeSz sources were investigated us¬ 
ing an ISW map constructed based upon the photometric 
luminous red galaxy (LRG) sample at redshift = 0.5 
(Granett, Neyrinck & Szapudi 2009). Although the map 
contained hot and cold spots corresponding to the super¬ 
structures, the amplitude of the effect was slight, with only 
^rj^ciuster _ ^rj^void ^ Q Qg ^ 0.1/iK. This result seemingly 

contradicted the claim that the lO/xK anisotropy was arising 
from the linear ISW effect. To resolve this question we are 
motivated to improve the fidelity of the ISW reconstruction 
by building it using a galaxy sample with measured spectro¬ 
scopic redshifts. 

We are now in a position to carry out these investiga¬ 
tions using the SDSS-HI Baryon Acoustic Oscillation Sur¬ 
vey GMASS spectroscopic sample which covers the sky area 
and redshift range of the original GraNeSz analysis (Alam 
et al. 2015; Bolton et al. 2012). In general, the BOSS sam¬ 
ple is interesting for ISW studies because, due to the volume 
dependence, the sensitivity to ISW is expected to peak at 
redshift z ~ 0.5 (Afshordi 2004). Analyses have shown the 
cross-correlation measurements of the luminous red galaxy 
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Figure 1. Redshift distribution of the CM ASS samples In the 
North and South Galactic cap fields. The vertical dashed lines 
indicate the redshift ranges used in this study. 

sample in the North Galactic cap show an excess signal (Ho 
et al. 2008; Giannantonio et al. 2008; Granett, Neyrinck & 
Szapudi 2009), despite low significance (Sawangwit et al. 
2010). The complete photometric BOSS sample was anal¬ 
ysed by Hernandez-Monteagudo et al. (2014) and Giannan¬ 
tonio et al. (2014). 

Galaxies in the GMASS sample have spectroscopic red¬ 
shift measurements, as opposed to photometric redshifts 
in the photometric LRG samples used previously. While 
precise redshift measurements are not necessary for cross¬ 
correlation measurements using the projected galaxy density 
field, they are necessary to reconstruct the three-dimensional 
potential. 

We begin this article by detailing the spectroscopic 
GMASS sample in Sec. 2 including our corrections for the 
known systematic inhomogeneities. In Sec. 3 we then de¬ 
scribe the procedure to estimate the underlying density field, 
the three-dimensional gravitational potential and the linear 
ISW map. We validate the ISW map using the template fit 
statistic against the Planck GMB temperature map. In Sec. 
4 we carry out the consistency check between the ISW map 
and the angular galaxy-GMB cross correlation function. In 
Sec. 5 we characterise the GraNeSz super-structures by fit¬ 
ting the three-dimensional density profiles. In Sec. 6 we ad¬ 
dress our main concerns, ISW sources. We characterise the 
three-dimensional properties of the GraNeSz structures us¬ 
ing the GMASS sample and estimate their contribution to 
the linear ISW effect using the ISW map. For comparison, 
we consider the publicly available void catalogue by Sutter 
et al. (2014). We also construct a new super-structure cata¬ 
logue based upon the local peaks and troughs of the three- 
dimensional gravitational potential that should be ideal ISW 
sources. We conclude in Sec. 7 with suggestions for future 
investigations. 

2 DATA 

2.1 Galaxy sample 

We use the constant stellar mass (GMASS) galaxy sample 
from the Sloan Digital Sky Survey-HI (SDSS-HI) Baryon 
Oscillation Spectroscopic Survey (BOSS) Data Release 12 
(Alam et al. 2015; Bolton et al. 2012). We extract the 
photometric parent catalog from the targeting catalog for 
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Figure 2. The mean number density of the CMASS sample with 
respect to a uniform random catalogue as a function of stellar 
density (left), seeing (center) and Galactic extinction (right). The 
dashed line shows the systematic trends present in the catalogue. 
After applying corrective weights we find the trends indicated by 
the solid lines. The shaded regions indicate the 68% uncertainty 
range for Poisson counts. 



Figure 3. The CMASS completeness maps for the North and 
South galactic caps. The top frames show the systematic com¬ 
pleteness correlated with stellar density, seeing, extinction and 
galaxy surface brightness. The bottom frames show the spectro¬ 
scopic completeness. 


the ninth and final run. This provides a uniform photo¬ 
metric selection over the full survey. The targets are se¬ 
lected with the CMASS target selection flag and fibre mag¬ 
nitude limit (boss_targetl&B0SS_CMASS>0). We then look 
up the redshifts of these sources in the spectroscopic objects 
database selecting all primary objects (scienceprimary=l). 
We match the photometric and spectroscopic catalogs using 
position on the sky with a match radius of 1 arcsec. We take 
the BOSS galaxy redshift estimate z_noqso when success¬ 
fully measured (zwarningjnoqso=0). If the BOSS redshift es¬ 
timate is not available, we take the legacy redshift measure¬ 
ment as long as it was successfully measured (zwarning=0). 

The CMASS sample covers two discontiguous fields in 
the Northern and Southern Galactic caps. We use the sur¬ 
vey mask that characterises the CMASS targeting catalogue. 
We extract the polygons from the mask that correspond to 
the BOSS survey runs (2,5,7,8 and 9) based upon the chunk 
identifiers. We exclude sources that fall within the photo¬ 
metric quality, bright star and centre post masks. We use 
the Mangle code (Swanson et al. 2008) to carry out these 
tasks. 


The spectroscopic redshift distribution of the final cata¬ 
log is shown in Fig. 1. When estimating the underlying den¬ 
sity field we combine two subsamples with redshift ranges 
(I) 0.43 < 2 : < 0.55 and (II) 0.55 < 2 : < 0.65 described in 
Table. 1. We consider separately the Northern and Southern 
Galactic cap fields. 


2.2 Survey completeness 

The parent catalog from which CMASS targets were selected 
is known to be inhomogeneous. Seeing conditions, stellar 
density and Galactic extinction were found to impact the 
number density of galaxies with a dependence on galaxy 
surface brightness (Ho et al. 2012; Ross et al. 2012). Here, 
we quantify those effects and compute corrective weights. 

We first estimate the spectroscopic sampling rate. The 
number of sources in the photometric target sample is given 
by Ntarg of which Nobs have measured redshifts. After clas¬ 
sifying as stars or galaxies we have Nobs = ^gai + Nstar- 
We then define the spectroscopic completeness for galaxies 
as Cspec = Ngal/{Ntarg — Nstar)- Thc complctcness is Com¬ 
puted for each galaxy by counting the number of sources 
within a radius of R = 30 arcmin. From the completeness 
we derive a weight for each galaxy: Wspec,i = ^/cspec,i- The 
map of the spectroscopic completeness is shown in Fig. 3. 
At redshift 2 : = 0.5, 30 arcmin on the sky corresponds to 
a comoving separation of ll/i“^Mpc, so we do not expect 
effects correlated on the scale of the detector to influence 
large-scale modes in the density field. 

To address possible large-scale systematic effects, we 
construct maps of each potential source. We selected stars 
as in Anderson et al. (2014) by taking sources with 17.5 < 
Mered < 19.9 from the SDSS star database table to create an 
angular density healpix map (Gorski et al. 2005). We build a 
map of the seeing taking the z-band FWHM of the PSF from 
the target photometric catalog averaging the seeing values 
within Healpix cells. We create a galactic extinction map in 
the same way using the i band extinction. We also use the 
2” fibre magnitude as a proxy for galaxy surface brightness. 

We compare the density of galaxies to the density of 
unclustered random points distributed over the survey area. 
Each random point has an associated stellar density, seeing 
and extinction assigned based upon the sky coordinate and 
a surface brightness drawn independently from the observed 
distribution. Fig. 2 shows the normalised galaxy density with 
respect to randoms as a function of each potential system¬ 
atic. There is a deficit of sources in regions with high stellar 
density or poor seeing while we find an excess of sources in 
regions with high values of reddening. 

The completeness depends on a four-dimensional pa¬ 
rameter space of fibre magnitude, stellar density, seeing and 
extinction. We develop a non-parametric approach using the 
K-means algorithm to adaptively bin the data. We use the 
algorithm to construct K = 500 adaptive bins each asso¬ 
ciated with roughly 2000 unique galaxies. This is carried 
out using a rank-order distance metric to standardise the 
parameters. 

We compute the completeness in each bin by counting 
the weighted number of galaxies with respect to the number 
of randoms: Csys = 'w;spec,i<^/W- The completeness is 

then assigned to each galaxy in the bin. The map of the 
systematic completeness is shown in the top row of Fig. 3. 
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Sample 


Number 

Median z 

Galaxy bias 

0.43 < z < 0.55 

North 

309900 

0.503 

1.94 ±0.01 


South 

109793 

0.499 

1.97 ±0.01 

0.55 < z < 0.65 

North 

218311 

0.590 

2.15 ±0.01 


South 

77872 

0.590 

2.18 ±0.02 


Table 1. Galaxy samples and best-fit bias parameters. The bias 
parameters are derived after fixing the amplitude of the matter 
power spectrum with erg = 0.8 under the fiducial cosmology given 
in the text. 


After weighting the density field by the inverse of the 
systematic completeness, we find the trends given by the 
solid red curve in Fig. 2. The correlations are reduced al¬ 
though we still see a significant dependence at the parame¬ 
ter extremes. However, the last bin in the plots is an average 
over 150 cells, so it represents only 1% of the total area. The 
trends shown in Fig. 2 are in good agreement with those 
published in Ross et al. (2012) which gives us confidence in 
our independent methodology. The total selection function 
is given by the product of the spectroscopic and systematic 
completeness fractions: Ci = CspecG x Csys,i- 


2.3 CMB maps 

The second data release from the Planck mission provides 
both temperature and polarisation maps (Planck Collabo¬ 
ration et al. 2015a). We use the foreground cleaned SMICA 
temperature map downgraded to 1° resolution (rigide = 64) 
(Planck Collaboration et al. 2015b). The polarisation data 
may also be used to help separate the late-time and pri¬ 
mordial temperature anisotropies. However, as of the latest 
release, the polarisation data is not available on large scales, 
harmonic ^ < 40, making the gain in the joint analysis negli¬ 
gible (Planck Collaboration et al. 2015c). Thus, we proceed 
with the analysis using only the temperature maps. 


3 ISW MAPS 

A prediction of the linear ISW signal may be obtained from 
the galaxy density field by employing the Poisson equa¬ 
tion (Granett, Neyrinck & Szapudi 2009). We use a Wiener 
filter to estimate the underlying matter density field and 
correct for survey completeness. The Wiener filter gives 
the minimum-variance estimate of the density field. In the 
Bayesian formulation it corresponds to the maximum a- 
posteriori solution assuming a Gaussian likelihood and prior, 
assumptions that hold for the galaxy field on large scales. 

Our method is based on the Gibbs sampler developed 
to reconstruct the underlying density field of the VIPERS 
galaxy sample (Granett, B. R. et al. 2015) with adjust¬ 
ments to accommodate the CMASS sample. We upweight 
the galaxies using the inverse of the completeness as de¬ 
scribed above. We divide the CMASS sample into two sub¬ 
samples by redshift: 0.43 < z < 0.55 and 0.55 < z < 0.65. 
The subsamples are combined to jointly estimate a single 
density field, but we allow the relative galaxy biases to vary. 
The weighted number of galaxies from each subsample is 
assigned over a comoving cartesian grid with cubic cells 


10/iMpc“^ on a side. We process separately the Northern 
and Southern Galactic caps. The North is embedded into a 
cube with 256 cells on a side while the South fits in a slightly 
smaller box with size 200 cells. The number of galaxies in 
a given cell indexed by i may be related to an underlying 
over-density field 

Ni^WiNi{l + 5g^i). ( 1 ) 

The mean number density of the cell is Ni and the selection 
function is wi. The selection function is estimated from the 
normalised number of random points falling in the given 
cell. It includes the angular mask. The redshift distribution 
is contained in Ni. The expected shot noise in the cell is 
erf = awiNi, where a accounts for the completeness weights. 

As in Granett, B. R. et al. (2015) our data model in¬ 
cludes a constant galaxy bias and the redshift distribution 
is left free. However, in this analysis we fix the power spec¬ 
trum to the fiducial model computed with GLASS to reduce 
the convergence time of the Gibbs sampler. The nonlin¬ 
ear power spectrum is computed using Halofit (Takahashi 
et al. 2012) with the fiducial parameter set {h = 0.6774, 
Qm = 0.3089, cTs = 0.80) at effective redshift Zref = 0.55). 
Use of the non-linear power spectrum aids in accurately 
recovering small-scale structures. Instead if we used the 
purely linear power spectrum the field would be unneces¬ 
sarily smoothed. We apply the Wiener filter isotropically in 
redshift space and so we use the monopole of the power spec¬ 
trum estimated using the dispersion model with distortion 
parameters 13 = 0.375 and av = 400km/s (Peacock V Dodds 
1994; Gabre V Gaztanaga 2009). 

We find that the Gibbs sampler converges after 10 steps. 
To produce realisations of the density field we ran 7 chains, 
each with 500 steps. The constraints on the bias parameters 
are listed in Table 1. 

From the density field we may estimate the gravita¬ 
tional potential using the Poisson equation in Fourier space, 

e^k) = -^H^nm{l + z)6ik). ( 2 ) 

In doing so, we scale S with the linear growth factor to a 
common epoch: S(zref) = 3{z)D{zref)/D{z). This compu¬ 
tation is carried out via the FFT algorithm. 

The ISW signal is an integral over the time-derivative 
of the potential along the line-of-sight, 

Tisw = -Tcmb-^ J (3) 

However, this is simplified using the linear prediction, that 
is, d^jdr oc d{l + z)D{z)/dr^ to obtain, 

Tisw(0,(j)) — —Tcmb^jy 
c G 

i 

which represents a projection over the potential field with a 
given kernel. After projecting the potential along the line- 
of-sight of the North and South samples we find the ISW 
maps shown in Fig. 4. In Appendix A we compare this map 
with the reconstruction from Granett, Neyrinck V Szapudi 
(2009) which was based upon a photometric luminous red 
galaxy sample over a similar redshift range in the Northern 
field. 


H“ Zre f)D{Zref) 
\d{l+z)D{z)] 
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Figure 4. Top row: CMB maps after low-pass filtering i < 20. Bottom row: ISW maps. Cross symbols mark the angular posistion of 
local maxima of the gravitational potential and open circles mark local minima more than 10° from survey boundaries. The locations of 
GraNeSz supervoids and superclusters are indicated by the large open circles and crosses. The CMASS survey boundary is overplotted 
in black over the CMB temperature maps. 


3.1 Template fit 

The observed CMB temperature map may be written as 
a sum of components including the primary temperature 
Tprimary^ thc ISW signal Tjsw and the uncorrelated fore¬ 
ground contribution Tforeground 


Tcmb — Tprimary T Tforeground + XTisw (5) 

The maximum likelihood estimator for the amplitude 
parameter A is 


c TcmbC ^Tisw 

~ TiswC-^Tisw ■ 


( 6 ) 

The variance is given by 



= {TiswC ^Tisw) 


( 7 ) 

The covariance matrix C is given by C = 

{(Tprimary ~b 


Tforeground)‘^)- Wc usc the bcst httiug tcmpcraturc power 
spectrum model from Planck computed using CLASS. In 
pixel space, the covariance matrix is 

2/ _i_ 1 

= 12 -^CiW^Piicoseij), ( 8 ) 

where Pi is the Legendre polynomial and Oij is the angular 
separation between two healpix cells indexed by i and j. 
The window function Wi is set by the Healpix pixelization 
at riside = 64. At this resolution we neglect the beam and 
instrument noise. 

A correction to the covariance must be added to ad¬ 
dress the fact that we estimate the mean of the distri¬ 
bution from the sample itself. The covariance is = 


{{xi — x){xj — x)) where x is the mean pixel value within 
the survey: x = J]). Xi. Expanding, this becomes = 

(xiXj) - Y.p{XaXp)- 

The results of the template ht are presented in Table 2. 
Considering the Northern map, we hnd that the ISW ampli¬ 
tude A = 0.41 with signal-to-noise 2.1. The Southern map 
gives no indication of a signal with an amplitude consistent 
with zero. In combination, we hnd A = 0.33 zb 0.17 giving a 
one-tailed p-value p = 0.03. 


3.2 Discussion 

The template ht amplitude gives the degree of agreement 
between the huctuations observed on the CMB and those 
expected from the galaxy held through the linear ISW effect. 
If the ingredients of the model are correct we expect an 
amplitude of unity. The value less than 1 means that we 
have over-predicted the huctuations on the CMB. Physically 
we may relate this to the amplitude of galaxy huctuations 
through the bias parameter implying a galaxy bias larger 
than anticipated. The amplitude A = 0.4 may be reconciled 
by increasing the bias by a factor of 2.5. 

Previously, Granett, Neyrinck & Szapudi (2009) found 
a template ht amplitude greater than 1. In that work, it was 
found that the amplitude depended on the number density, 
redshift distribution and redshift uncertainty of the galaxy 
tracers and furthermore, that shot noise biases the ampli¬ 
tude low. Thus at this point, we do not draw conclusions 
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ATisw Template fit 
Field A cr A/cr 

1-tail p 

North 

0.41 

0.20 

2.1 

0.02 

South 

0.09 

0.34 

0.3 

0.4 

Combined 

0.33 

0.17 

1.9 

0.03 


Power spectrum fit 


Field 

A 

(7 

A/cr 

1-tail p 

North 

2.5 

1.1 

2.2 

0.01 

South 

2.1 

1.9 

1.1 

0.1 

Combined 

2.4 

1.0 

2.4 

0.01 

Clean North 

0.7 

1.1 

0.6 

0.3 

Clean South 

1.8 

1.9 

0.9 

0.2 

Clean Comb. 

0.9 

1.0 

0.9 

0.2 


Table 2. The best-fit model template amplitudes for the ISW 
map and angular power spectrum measurements. 


from the value of the amplitude and take the signal-to-noise 
ratio as confirmation of the validity of the map. 


4 ANGULAR CROSS-CORRELATION 

A standard probe of the ISW effect is the cross-correlation 
of galaxy density fields and CMB temperature maps. For 
SDSS luminous red galaxy data sets, previous studies found 
a moderately significant 2 — 2.5cr correlation with ampli¬ 
tudes higher than expected in the ACDM model (Ho et al. 
2008; Giannantonio et al. 2008; Granett, Neyrinck Sz Sza- 
pudi 2009). Other studies reported ISW signals consistent 
with ACDM (Hernandez-Monteagudo et al. 2014; Giannan¬ 
tonio et al. 2014). 

We select GMASS galaxies in the redshift range 0.43 < 
z < 0.65 and compute the projected density field on a 
Healpix grid with resolution 1° (riside = 64). The selection 
function is estimated by assigning particles from the random 
catalogue in the same fashion. We then normalise the galaxy 
count map by the density randoms in each cell, setting to 
zero cells that are less than 50% complete. 


4.1 Auto power-spectra 

We compute harmonic-space power spectra, = 

Wirn^Ym)^ using the SpIGE code (Ghon et al. 2004) which 
properly accounts for the survey geometry. We first com¬ 
pute the auto-power spectrum of this map as a function of 
spherical harmonic index shown in Fig. 5. We mea¬ 

sure consistent values in our bins of width Ai = 10 in the 
Northern and Southern hemispheres. We thus analyse mul¬ 
tipoles ^ < 80 since the theoretical S/N for an ISW signal 
in the ACDM model already converges at ~ 1.7 for £ ^ 60 
(Hernandez-Monteagudo et al. 2014). For comparison, we 
overplot a model computed under best-fit Planck cosmo¬ 
logical parameters {h = 0.6774, = 0.3089, crs = 0.80) 

at effective redshift Zref = 0.55. The constant galaxy bias 
bg = 2 provides a good fit. The systematic weights alleviate 
the excess power on large-scales. 

The fiducial model is given by the projection of the 


rr {h ^Mpc) 

Supervoids 
a q 

aios (h ^Mpc) 

35.9 ±3 

1.7 ±0.15 

2.6 ±0.4 

76 ±9 

rs {h~^Mpc) 

Superclusters 

Sc q 

(7los (/i“^Mpc) 

9.0 ±3 


0.64lo;| 

79 ±9 


Table 3. Best fit parameters for GraNeSz supervoid and super¬ 
cluster mean profiles with 68% confidence limits. 


power spectrum in the plane-parallel limit: 


cr = 


J dr r'^n{r)^ bg{r) 


D\r) 

DRZref) 


I + 1/2 


Zref 


(9) 


The redshift distribution is normalised such that 
J dr r^n(r) = 1 , b(r) is the galaxy bias, D(r)/D(zref) is 
the growth factor that linearly evolves the field to redshift 
Zref and P is the non-linear power spectrum estimated at 
redshift Zref- We assume that the galaxy bias is constant 
over the redshift range of interest. 


4.2 Cross power-spectra 


Similarly, the expected angular power spectrum in the flat- 
sky limit, see Afshordi (2004), is an integral over the line- 
of-sight distance r. 


cr 


Tcmb 





(/ + 1 / 2)2 
d(l + z)Di{z) 
dr 


P 




( 10 ) 


where, n{r) = is the galaxy selection function 

normalized by f r^ growth factor and 

P (^k = is the matter power spectrum at redshift 

Zref • 

In Fig. 5 we show the cross correlation between the 
galaxy density field and the GMB. We allow the amplitude 
of the model to be free with parameter A and find the best-fit 
value and uncertainty. We introduce the amplitude parame¬ 
ter A which has a distinct interpretation from the amplitude 
A used in the template map fit. The value A = 1 would im¬ 
ply consistency with the AGDM expectation with the best- 
fit galaxy bias and A > 1 can be reconciled by increasing 
the galaxy bias. The fit results are summarised in Table 2. 
While systematic differences between the two hemispheres 
may affect our measurement (Giannantonio et al. 2014), we 
find consistent ISW signals for the Northern, Southern, and 
combined maps. 

As a consistency test of the ISW map derived in Sec. 
3, we subtract it from the Planck GMB temperature map 
and refer to this as the ISW cleaned map. Repeating the 
cross-correlation measurement, we find that the amplitude 
drops to A = 0.9zbl.0 when the galaxy density map is 
cross correlated with the ISW cleaned GMB map. The match 
between the two methods is better for the Northern field, 
while a residual remains in the South at the Icr level. 
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Figure 5. Left: the angular auto-power spectra of BOSS galaxies. The markers with error bars give the mean value in a bin of width 
Ai = 10. We show the measurement with and without the corrective weights for foreground systematics. Over-plotted with the red solid 
curve is a fiducial power spectrum computed with CLASS and halofit. We find that the data is consistent with a bias bg = 2.0. Right, the 
cross galaxy-CMB temperature power spectrum. The fiducial model is over-plotted and we find a best fitting amplitude A = 2.4 ± 1.0. 
After subtracting the ISW map from the CMB, the cross-correlation drops, shown by the arrows. 


5 CHARACTERISATION OF GRANESZ 
SUPER-STRUCTURES 

5.1 Supervoid and supercluster profiles 

The spectroscopic CM ASS sample allows us to measure the 
shape and sizes of GraNeSz super-structures which were 
originally identified in a photometric redshift sample. We 
begin by measuring the number density of CM ASS galax¬ 
ies around the positions of GraNeSz super-structures. The 
catalogue includes 50 supervoids and 50 super clusters. We 
measure the mean density profile by taking the ratio of the 
number of galaxies to the normalised number of unclustered 
random points in a given separation bin. The number of 
points in the random catalogue is 17 times the total number 
of gsalaxies. Fig. 6 and 7 show the mean number density 
of the 50 supervoids and 50 superclusters as a function of 
distance parallel and perpendicular to the line-of-sight. We 
find that the profiles are not localised in the line-of-sight 
direction. This may be due to error in the individual cen¬ 
tres estimated with photometric redshifts. To interpret these 
plots we will fit simple void and cluster profile models that 
includes uncertainty in the structure centre. 

We estimate the variance of the profile measurements by 
generating many realisations centred on random positions in 
the survey. First, considering either the supervoid or super¬ 
cluster catalogue, 50 random points are drawn from a ran¬ 
dom catalogue over the Northern field survey volume. Red- 
shifts are assigned from the true redshifts in the superstruc¬ 
ture catalogue. Then, the mean profile is measured from 
these positions. We used 1000 realisations of the process to 
estimate the diagonal elements of the covariance matrix for 
both supervoid and supercluster catalogues. 


5.2 Void and cluster models 

We use the following analytic form for the void profile with 
scale radius Vy and shape parameter a (Hawken et al, in 



Figure 8. The GraNeSz supervoid profile (left) projected parallel 
to the line-of-sight, (center) projected perpendicularly to the line- 
of-sight, (right) the deconvolved three-dimensional model. The 
measurements are compared with the best-fit model (dark solid 
curve) and samples taken from the Markov chains (light blue 
curves). 




Figure 9. The same as Fig. 8, but for GraNeSz superclusters. 


prep.): 




exp 



( 11 ) 


The void is permitted to be elongated along the line-of-sight 
with ellipticity parameter g, and we define: 

= rlerp+rlar/q'^- ( 12 ) 


Superclusters are known to have irregular morphologies 
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Perpendicular separation (h 'Mpc) 


Figure 6. Mean GraNeSz supervoid profile: (a) the normalised number density of CMASS galaxies measured in bins perpendicular and 
parallel to the line-of-sight (bin size is 5/i“^Mpc), (b) the best-fit model density profile of the void, (c) the model convolved with the 
best-fit Gaussian line-of-sight position error showing how it would appear in observations, (d) the residual in units of standard deviation 
after subtracting the convolved model from the data. 



Perpendicular separation (/? ^Mpc) 


Figure 7. Same as Fig. 6, but for the GraNeSz supercluster catalogue. 


(Einasto et al. 2014) and may not be well approximated 
by spherical symmetry. Nevertheless we apply the Navarro, 
Frenk & White (1996, NEW) analytic form to characterise 
the mean supercluster profile: 


P(r) ^ ^ G 

P f{l + f/rsf 


(13) 


The supercluster profile is also allowed to be elongated along 
the line-of-sight with parameter q. 

We allow for a scatter in the line-of-sight position by 
introducing a convolution with a Gaussian kernel charac¬ 
terised with scale aios’- 


pzG. 


perpf I par 


) = 


1 



(y-rpar)^ 

los pij'perp'j y)dy. 

(14) 


We ht the model parameters to the measurements using 
a Markov chain monte carlo method implemented with the 


Emcee package^ (Foreman-Mackey et al. 2013). The ht is 
carried out using the diagonal covariance matrix constructed 
as described above. 

We give the medians of the parameter chains along with 
the 68% conhdence levels in Table 3. In Figs. 6 and 7 we 
show the mean anisotropic prohles of the voids and clus¬ 
ters along with the best-ht models with their residuals. The 
ht is carried out on the model convolved with the Gaus¬ 
sian kernel (Eq. 14) but we also show the prohle before the 
convolution step to represents the underlying galaxy distri¬ 
bution. Figs. 8 and 9 compare the model with perpendic¬ 
ular and parallel projections of the prohles along with the 
model three-dimensional radial prohle. The goodness of ht 
is xloid = 0-93 and x%uster = 0.99 per degree of freedom for 
the supervoid and supercluster models respectively. 


^ http://dan.iel.fm/emcee/ 
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Radial distance {h ^Mpc) 

Figure 10. Radial density profiles of BOSS super-structures. Top 
row: local minima and maxima of the gravitational potential. 
Middle row: GraNeSz supervoids and superclusters. The best- 
fit models for the three-dimensional profile parallel (long-dashed 
curve) and perpendicular (dot-dashed curve) to the line-of-sight 
are shown. Bottom row: supervoids from the Sutter DRIO anal¬ 
ysis. Plotted are profiles in spherical shells as well as cylindrical 
shells along the line of sight within r = 200/i“^Mpc to show a 
representation of the structure in photometric redshift space. 

In principle we could recenter each structure along the 
line-of-sight to correctly superpose the profiles. However, we 
found that due to the sparseness of the spectroscopic cat¬ 
alogue, the true centres could not be accurately estimated 
and so the uncertainty in the position can only be corrected 
for statistically on the mean profile. 


5.3 Discussion 

The model fits demonstrate that the supervoid and super¬ 
cluster samples both represent large structures. The voids 
reach to R = 50/i“^Mpc and the scale radius for superclus¬ 
ters is Vs = 9/i“^Mpc. 

The mean supercluster profile is much broader than 
that which would be expected for a typical massive clus¬ 
ter. We find that the radius that encloses a mean density 
200pcrit is r 2 oo = l.l/i“^Mpc giving a concentration index 
C 200 = rs/r 2 oo of C 200 = 0.12 while massive clusters have 
concentration C200 ~ 3 (Umetsu et al. 2015). Given the low 
concentration index, the mass within r 2 oo is modest: we find 
M 200 = 1.5xlO^^/i“^M0.Ifwe instead integrate to the scale 
radius Ts we find M{r < r^) = 4.7 x lO^^/i^^M©. For a typi¬ 
cal massive cluster these masses are M200 = 1 x lO^^/i^^M© 
and M(r < Vg) = 2x lO^^/i^^M©. The scale of the superclus¬ 


ters suggests that they are not virialized structures. Despite 
this, we find that the NFW profile characterises the mean 
profile well. 

The best-fit void model is elongated along the line-of- 
sight with axis ratio q = 2.6 zb 0.4. We found that a spher¬ 
ically symmetric model is unable to fit the data even with 
the free line-of-sight smoothing parameter. We explored this 
peculiarity with more general void profile forms, but found 
that even a spherically symmetric step function profile is un¬ 
able to reproduce the measurement. On the other hand, the 
supercluster model is not strongly elongated: we find only a 
marginal squashing with q discrepant from unity at the Icr 
level. 

How can we interpret the apparent elongation of the 
supervoids? This is likely to be a real effect arising from 
the void-finding algorithm. GraNeSz identified the super¬ 
voids and superclusters in a photometric redshift catalogue 
with redshift error cr^ ~ 0.05 ~ 100/i“^Mpc. We suggest 
that when applied to a photometric sample the algorithm is 
most sensitive to under-dense tubes that are pointed along 
the line-of-sight. The effect is important for voids because 
they are large, irregularly shaped structures. Superclusters 
are also known to have wall-like or filamentary morphologies 
(Einasto et al. 2014); however, we find that the structures 
in the supercluster catalogue do not have a preferred align¬ 
ment. 

The void and cluster models incorporate the posi¬ 
tion uncertainty through the smoothing parameter aios- We 
found that the parameter is required to fit both the super¬ 
void and supercluster measurements; moreover the scales 
determined from the two models are in excellent agreement. 
Interpreting this effect as arising from redshift error, the 
scale aios = 76/i“^Mpc corresponds to a redshift error of 
az = 0.034. 


6 COMPARISON OF BOSS 
SUPER-STRUCTURES 

6.1 Local extrema of the gravitational potential 

We construct a new superstructure catalogue defined by 
the positions of local extrema of the gravitational poten¬ 
tial field. The starting point is the Wiener estimate of the 
gravitational potential derived through the Poisson equa¬ 
tion, Eq. 2. We perform the computation on a grid with cell 
size 10/i“^Mpc, thus we may resolve structures on scales 
R > 20/i“^Mpc. To find the largest structures, we convolve 
the field with a three-dimensional Gaussian filter with full- 
width half-maximum 60/i“^Mpc. The local extrema of the 
field then correspond to cells that are more extreme, either 
greater or less than, their six adjacent cells. The positions 
of the extrema are discretised by the grid. This introduces 
an uncertainty of 10/i“^Mpc which corresponds to an angu¬ 
lar scale of ~ 30 arcmin; however, this uncertainty in the 
structure positions is much less than the smoothing scale. 
Minima of the potential correspond to over-densities, while 
maxima correspond to under-densities in the galaxy field. 
We exclude objects with centres less than 10° from the sur¬ 
vey boundary. We find 136 and 25 supervoids in the North 
and South, respectively, and 166 and 49 superclusters in the 
North and South, marked in Eig. 4. 
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6.2 Radial density profiles 

In Fig. 10 we compare the radial density profiles of three 
super-structure catalogues: GraNeSz, peaks and troughs of 
the potential and the DRIO spectroscopic void catalogue by 
Sutter et ah (2014). The projected profiles are computed 
by averaging in cylinders along the line-of-sight within a 
sphere of r = 200/i“^Mpc. In the case of the GraNeSz struc¬ 
tures the void centres are given by the photometric redshift 
estimates, but we overplot the deconvolved best-fit model 
found in Sec. 5.2. The Sutter voids were selected in the red- 
shift range 0.5 < z < 0.6 and have radii R > 50/i“^Mpc. 
We found that for voids with radius less than 40/i“^Mpc 
the structures in projection are overdensities, a fact pointed 
out by (Nadathur & Hotchkiss 2014; Gai et al. 2014). We 
find that the GraNeSz voids are deeper in projection than 
the Sutter voids, reinforcing the fact that photometric and 
spectroscopic void finders can be sensitive to different types 
of structures. 

The superstructures identified from local extrema of 
the gravitational potential appear significantly larger in ex¬ 
tent. The supervoids appear significantly underdense, deeper 
than GraNeSz supervoids and match the depth of Sutter 
et al. (2014) catalogue, but lacking the compensated ring 
structure. Moreover, in projection they appear deeper than 
the GraNeSz voids. The superclusters are found to be sig¬ 
nificantly more dense than GraNeSz superclusters, both in 
three-dimensions and projection. 


6.3 Radial temperature profiles 

We now examine the imprint of these structures in our model 
ISW map and in the Planck temperature map. We will use 
the ISW map as a template to match to the GMB data. 
The radial ISW temperature profiles are shown in the left 
column of Fig. 11 and the temperature profiles on the GMB 
are shown on the right. The shaded regions show the 68 
and 95% confidence intervals. For the ISW map, the con¬ 
fidence interval was determined by measuring the profile 
around randomly distributed points within the survey mask. 
We use 1000 random realisations to compute the variance of 
the profile in this way. On the right, for the variance on the 
GMB temperature, we keep the coordinates of the structures 
fixed and generate mock GMB skies using the Healpix Syn- 
fast code. Here we also use 1000 realisations and compute 
the covariance matrices for the profiles. 

The amplitude of the profile is affected by long- 
wavelength fluctuations that effectively modulate the back¬ 
ground level. To remove this variance we apply a high-pass 
filter to both the ISW and GMB maps. We use a filter with 
a cosine cut-off at ^ = 10, corresponding to 20° scales. In 
presenting the profile, we further subtract the mean tem¬ 
perature at 10° scales. We do not expect such filtering to 
affect a well-localised signal arising from individual voids or 
clusters. 

We see that each super-structure sample is correlated 
with the ISW map to some degree. The peaks and troughs 
of the potential have the largest amplitude, which can be 
expected by their construction. The GraNeSz sources do not 
correspond to the most extreme features in the map, but 
show significant correlation. On the other hand, the Sutter 
et al. (2014) sample shows only a weak correlation with the 


Sample 

N 

A 

cr 

Xja 

GraNeSz supervoids 

50 

6.30 

2.59 

2.4 

GraNeSz superclusters 

50 

4.51 

3.57 

1.3 

Potential minima 

215 

0.23 

0.59 

0.4 

Potential maxima 

161 

0.61 

1.05 

0.6 

Sutter DRIO 

112 

1.23 

5.39 

0.2 


Table 4. ISW source template fit significances. 
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Figure 11. Mean temperature profiles of super-structure cata¬ 
logues. Left: the profile of structures projected on the ISW map; 
right: the measured temperature profile on the CMB. The shaded 
regions show 68% and 95% confidence intervals. 


ISW map, which may be anticipated by their shallow profiles 
in projection. 

We use the ISW profile measured on the ISW map as a 
template to fit the GMB temperature profile. The template 
fit statistic and variance are computed as given by Eqs. 6 and 
7. The resulting significances are listed in Table 4. We find 
marginal detections of the GraNeSz supervoids and super¬ 
structures. The supervoid sample gives higher significance 
than the supercluster because there is better agreement be¬ 
tween the measured profile and the template shape. For the 
other samples the fits are consistent with no correlation. 


7 CONCLUSIONS 

The BOSS survey volume and sampling rate allow a de¬ 
tailed investigation of ISW sources at redshift z = 0.5. We 
have used the spectroscopic sample to further investigate 
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the GraNeSz super-structure catalogue and carried out a 
preliminary analysis of public catalogues. Our main findings 
may be summarised as follows. 

• We constructed a map of the ISW anisotropy sourced 
by CMASS galaxies and confirmed that the signal is present 
in the CMB temperature map at the 97% confidence level. 
The detection is not significant in the Southern field. This 
may be due to the small area and intrinsic variance of the 
signal or contamination by strong foreground systematics in 
the South. 

• We demonstrated that the signal is consistent with 
the measured galaxy density - temperature cross-correlation 
function with similar North-South differences. 

• We carry out a detailed characterisation of the GraNeSz 
super-structures. We confirm that they are indeed vacuous 
supervoids and extensive supercluseters, but they are not 
well localised in the line-of-sight direction. We model this 
uncertainty as a convolution of the profile. Fitting the model 
to the observed mean profile, we find that the supervoids are 
intrinsically elongated along the line-of-sight. We suggest 
that this feature arises from the void-finding algorithm run 
on a photometric redshift sample. 

• We construct a new super-structure catalogue based 
upon peaks and troughs of the gravitational potential. These 
structures have a much larger extent that the GraNeSz and 
Sutter et al. (2014) objects, and should be expected to give 
a strong linear ISW contribution. 

• We compare the mean radial density profiles of struc¬ 
tures and find that Sutter et al. (2014) voids are relatively 
shallow in projection owing to the compensated shell. The 
voids found with photometric redshifts do not show this 
morphology. 

• We construct ISW signal templates for the superstruc¬ 
ture catalogues based upon our ISW map. We use the tem¬ 
plate as a filter on the GMB to measure the ISW con¬ 
tribution, largely free of ‘look elsewhere’ bias. With this 
method we find marginal detections of the GraNeSz super¬ 
structures, but the other super-structure samples are con¬ 
sistent with no signal. 

We have resolved remaining questions regarding the na¬ 
ture of the GraNeSz super-structures by quantifying their 
morphologies; however, their apparent imprint on the GMB 
remains unexplained. The elongation of the voids may offer 
a clue: we expect that the extreme tail of the ISW signal is 
sourced by structures that happen to be elongated along the 
line-of-sight (Flender, Hotchkiss & Nadathur 2013). How¬ 
ever, the linear ISW signal that we estimate based upon 
the reconstructed potential does not account for the mea¬ 
sured temperature anisotropy. Moreover our analyses of 
other super-structure catalogues give results fully consistent 
with the expectation in the AGDM model. Thus we are un¬ 
able to confirm the existence of anomalous ISW sources. 

The approach we have developed here is fully empiri¬ 
cal: we can characterise the sources expected to contribute 
most strongly to the ISW signal. However there is also an 
important need for theoretical expectations that can be pro¬ 
vided by N-body simulations. The picture may come into fo¬ 
cus with future refinements of super-structure catalogues in 
BOSS. We also look to upcoming photometric and spectro¬ 
scopic surveys to improve the fidelity of the galaxy density 
field reconstructions on large scales. 


APPENDIX A: COMPARISON WITH 
PHOTOMETRIC ISW MAP 

We have probed the ISW signal over the same area of sky 
and in a similar redshift range as an earlier work by Granett, 
Neyrinck & Szapudi (2009). This previous study used a pho¬ 
tometric sample of luminous red galaxies in the redshift 
range 0.4 < Zphot < 0.6 based on photometric redshifts. 
Since the integrated linear signal depends on the redshift of 
the sources only weakly through the growth factor we expect 
that maps generated with photometric or spectroscopic red¬ 
shifts should be equivalent. The maps from the two works 
are shown in Fig. Al. By eye we can find some correspon¬ 
dence between the two maps for instance the hot spot in the 
top left at RA, dec. 210°, 50° is in common; however, the 
large-scale fluctuations do not appear to match up. 

For quantitative comparison we compute the correla¬ 
tion coefficient of the spherical harmonic decompositions of 
the two maps A and B indexed by t Vi = j 
plotted in the bottom panel of Fig. Al. On large scales we 
find an anti-correlation with r < 0. The correlation becomes 
positive at ^ > 10 (20° scales) with r ~ 0.6. The mean corre¬ 
lation between the two maps is r = 0.2. As a sanity-check we 
flip one of the maps in right ascension and repeat the mea¬ 
surement. This results in a reduced correlation coefficient as 
shown by the dashed curve in the figure. 

We find that the two maps are correlated on small scales 
at the 60% level but it is not perfect. This may be explained 
by the noise properties of the samples: the photometric sur¬ 
vey has a greater density of galaxies and so gives an estimate 
with reduced shot noise. Furthermore, the redshift distribu¬ 
tions of the samples do differ. 
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Figure Al. A comparison of the ISW map derived in this work 
(top frame) with the map estimated using a photometric sample 
of luminous red galaxies from Granett, Neyrinck & Szapudi (2009) 
(middle frame). The correlation coefficient between the two maps 
is plotted in the bottom frame as a function of harmonic multipole 
i (solid line). The dashed line shows the result after flipping one 
of the maps in right ascension. 
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